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Research  Summary 

The  research  performed  under  this  project  has  emphasized  an  improved 

characterization  of  the  spect .roscopic  properties  of  the  diatomic  mercury  halide 

molecules.  The  work  has  employed  mainly  conventional  methods  of 

high-resolution  emission  spectroscopy,  coupled  with  the  use  of  "cool”  Tesla 
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discharge  sources  containing  single  isotopomers  (e.g.  ~  Hg  Br)  .  The  use  of 
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isotopically  pure  ~  Hg  in  these  sources  has  led  to  considerable  simplification 
of  the  emission  spectra,  from  elimination  of  the  isotopic  blending  that  results 
from  vibrational  isotope  shifts  in  "natural"  HgX  sources.  The  primary  targets 
of  these  studies  have  been  the  B*X  emission  systems,  which  occur  in  the  visible 
region  of  the  spectrum  and  serve  as  the  lasants  in  the  mercury  halide 
"blue-green"  lasers.  However,  our  more  recent  work  has  included  reanalysis  of 
UY  emission  systems  in  these  molecules,  in  order  to  improve  the  description  of 
the  HgX  ground  states  at  low  v.  where  they  are  not  accessed  in  emission  from 
the  B  states  under  the  conditions  of  our  discharges  (and  of  HgX  lasers). 

These  studies  have  led  to  extensive  revisions  in  the  vibrational 
assignments  of  bands  which  occur  in  the  most  intense  region  of  the  HgX  B->X 
spectra,  which  are  the  regions  relevant  to  the  HgX  lasers.  One  outcome  of 
these  revisions  is  a  lowering  fby  about  lOfc)  of  the  estimated  dissociation 
energies  of  the  ground  ( X)  states  of  these  molecules.  Accordingly  the 
potential  curves  for  these  sta*'-  are  altered  from  their  previous  assessment, 
and  with  them  the  calculated  Franck-Condon  factors  for  B*X  emission  from  low  v' 
levels.  The  new  FCFs  are  more  in  keeping  with  the  observed  features  in  the 
laser  spectra,  although  there  is  still  no  simple  assignment  of  the  latter. 
Rather  they  appear  to  be  associated  with  chance  coincidences  of  many  rotational 
lines  in  the  densely  overlapped  spectra  of  the  multiple  isotopomers  occurring 
in  the  "natural”  HgX  molecules  employed  in  the  laser  experiments.  Several 
papers  have  been  published  on  this  work  (items  1,2, 6-8  under  "Publications" 
below)  and  more  are  anticipated  (see  Appendices  1  and  5). 

In  connection  with  our  study  of  the  HgX  emission  spectra,  we  have 
developed  new  methods  for  analyzing  spectra  of  heavy  diatomic  molecules.  In 
such  molecules  electronic  spectra  t'-nd  to  be  very  congested,  even  for 
relatively  low  discharge  temperatures.  This  congestion  is  a  direct  consequence 
of  the  small  vibrational  and  rotational  frequencies  which  always  occur  in 
diatomics  with  large  reduced  masses:  and  it  makes  a  rotational  analysis 
difficult  or  practically  impossible.  On  the  other  hand  the  vibrational 


structure  in  such  spectra  may  be  rich,  leading  to  quite  good  determination  of 
the  vibrational  constants.  We  hav<  found  that  it  is  possible  to  obtain  fairly 
reliable  potential  curves  for  electronic  states  involved  in  such  congested 
transitions,  bv  using  a  Morse  approximation  to  represent  the  unknown  left  or 
repulsive  branch  of  the  potential,  and  the  RKR  /  integrals  (which  require  only 
the  vibrational  constants)  to  calculate  the  width  of  the  potential  --  and  hence 
the  right  or  attractive  branch  —  as  a  function  of  v  and  the  vibrational  energy 
G  .  A  paper  describing  this  method  has  been  published  (item  3  under 
"Publications"  below). 

Another  method  we  have  developed  is  the  use  of  "mixed"  representations  f  r 
rotational  and  vibrational  constants  of  diatomics.  It  is  well  known  that 
polynomials  in  (v  +  1/2)  are  very  efficient  for  representing  B  and  G  da’ a  for 
diatomics  at  low-t o-moderate  v.  However,  near  dissociation  such  polynomials 
are  not  well  behaved  and  they  become  very  inefficient.  In  this  regime 
so-called  near-dissociation  expansions  (NDEs).  which  represent  the  energy  and 
rotational  constant  in  a  theoretically  sound  way  with  the  dissociation  licit  as 
reference,  become  very  efficient.  However.  NDEs  are  no  better  at  low  v  than 
are  polynomials  in  (v  +  1/2)  at  high  v.  Hence  we  have  experimented  with  miv -i 
representations  --  polynomials  for  low  v,  NDEs  for  high  —  and  have  found  th'-m 
to  perform  quite  well  in  our  test  calculations  on  HgBr(X) .  A  paper  describing 
this  method  has  been  published  (Publication  8)  and  represents  an  extension  of 
earlier  work  on  NDEs  (see  Publication  5). 

As  a  spin-off  of  our  work  on  efficient  representations  of  spectroscopic 
data  for  diatomics,  we  have  also  developed  new  methods  for  calculating  diatom’  • 
partition  functions.  These  methods  are  particularly  useful  for  dealing  with 
electronic  states  which  are  only  weakly  bound,  such  that  significant 
contributions  to  the  partition  function  (and  hence  to  the  thermodynamic 
properties)  come  from  levels  near  dissociation.  One  of  the  new  methods  can  l •• 
considered  exact;  it  employs  efficient  Gaussian  quadrature  schemes  to  evaluati 
the  semiclass i cal  eigenvalues  for  a  given  potential  curve,  and  then  ev.-uate- 
the  partition  function  as  a  direct  sum  of  the  Boltzmann  factors  for  the.-*- 
levels.  The  second  method  employs  classical  phase  integrals  to  evaluate 
rotation-vibration  partition  functions  about  two  orders  of  magnitude  faster 
than  the  semielassical  procedure,  with  a  reliability  of  about  lrt  when  kT  i> 
greater  than  about  twice  the  vibrational  frequency.  (This  condition  can  occur 
at  relatively  low  T  for  heavy  diatom  cs . )  These  methods  have  been  described  in 


two  published  papers  (Publi oat  ions  10  and  12.  included  as  Appendices  1  and  3  of 
this  report). 

In  an  extension  of  the  original  program  of  research  on  the  HgX  molecules. 
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we  have  also  worked  on  the  analysis  of  the  B+X  spectrum  of  uXe1JF.  This  work 
was  supported  by  DARPA  and  led  to  a  publication  of  an  updated  and  refined 
analysis  of  the  XeF  laser  spectrum  (Publication  11  and  Appendix  2). 

Continuing  Work 

We  have  accomplished  the  main  goals  of  this  project  —  namely  the 
reanalysis  of  the  emission  spectra  of  HgCl .  HgBr,  and  Hgl.  Brief  reports  on 
this  work  have  been  published,  and  the  full  papers  describing  the  work  on  Hgl 
and  HgBr  are  now  in  preparation  for  publication  in  the  Journal  of  Molecular 
Spectroscopy .  Most  of  the  tables  to  be  included  in  these  papers  have  been 
prepared  and  are  included  with  this  report  as  Appendices  4  and  5.  In  the  case 
of  HgCl  the  rotational  analysis  of  B-*X  remains  to  be  done.  This  work  should  be 
completed  by  September.  1985;  it  too  will  be  published  in  the  Journal  of 
Molecular  Spectroscopy. 

In  the  extension  work  on  XeF  we  have  also  succeeded  in  the  main  goal  of 
completing  the  rotational  analysis  of  the  B-*X  spectrum  of  XeIJF.  Here  too 
questions  remain,  particularly  as  regards  the  possibility  of  B-C  perturbations . 
which  should  occur  but  which  have  not  yet  been  identified  in  our  spectra. 
Again  a  brief  report  of  this  work  has  appeared,  and  the  full  publication  is 
planned  for  the  Journal  of  Molecular  Spectroscopy . 

Initially  we  intended  to  conduct  transient  absorption  experiments  on  the 
HgX  molecules  in  order  to  characterize  the  B~X  absorption  from  low  v' '  levels 
to  high  v '  levels  —  the  region  not  sampled  in  B+X  emission.  However  we  opted 
for  a  reanalysis  of  the  UV  emission  spectra  (primarily  C*X  and  D+X)  instead,  as 
these  transitions  nicely  determine  the  low-t>  regions  of  the  ground  states. 
This  work  has  led  to  greatly  improved  descriptions  of  the  ground  states, 
including  their  thermodynamic  properties  in  the  case  of  HgBr  (Publication  12). 
We  expect  to  publish  papers  on  the  partition  functions  and  thermodynamic 
properties  of  Hgl,  HgCl.  and  XeF,  also,  when  the  fundamental  spectroscopy 
papers  on  these  molecules  have  been  completed. 

There  is  one  element  of  our  project  in  which  we  have  failed  to  accomplish 
our  goals.  That  is  the  intended  study  of  pressure  broadening  of  rotational 
lines  in  the  HgX  and  XeF  B-*X  spectra.  To  date  we  have  been  unable  to  conduct 


very  high  resolution  measurements  with  our  Fabry-Perot  interferometer,  which  we 
had  hoped  to  couple  with  our  spectrometer  to  study  line  widths.  Although  we 
have  not  abandoned  this  goal,  it  is  clear  that,  success  in  this  area  will 
require  time  and  effort  beyond  the  realm  of  this  project,  so  this  work  will  not 
be  completed  in  the  foreseeable  future. 

In  another  area  —  the  study  of  pressure  and  temperature  dependence  of  the 
HgX  emission  spectra  —  we  have  done  some  work  (see  Publications  7  and  9)  but 
will  probably  not  complete  the  calculations  needed  to  estimate  vibrational 
relaxation  rates  and  the  /?-dependence  of  the  B-X  transition  moment  function. 
As  work  progressed  on  the  analysis  of  the  spectra,  it  became  clear  that  this 
aspect  of  the  study  could  provide  only  rough  estimates  at  best;  consequently  it 
was  relegated  to  low  priority. 

Personnel 

Three  graduate  students  —  K.  S.  Viswanathan.  J.  Gail  Ashmore,  and 
0.  Carlysle  Salter  —  have  been  employed  essentially  full  time  on  this  project . 
Dr.  Viswanathan  has  completed  the  requirements  for  his  Ph.  D.  degree  (thesis; 
"Part  1:  Spectroscopic  Studies  of  Charge  Transfer  Transitions  in  Iodine, 
Bromine,  and  Mercury  Iodide;  Part  2:  Nitrogen  Laser  Pumped  Dye  Laser,"  dune. 
1983).  and  has  taken  a  postdoctoral  appointment  in  the  Department  of  Chemistry 
at  the  University  of  Indiana.  Ashmore  and  Salter  are  now  in  their  fifth  vears 
and  expect  to  finish  their  theses  before  September,  1985.  Three  undergraduates 
have  also  been  associated  with  the  project:  Patrick  Berwanger.  Sue  A.  Davies, 
and  Stuart  D.  Henderson.  Berwanger  is  currently  employed  with  Exxon  in 
Houston,  Texas,  and  Davies  and  Henderson  are  in  their  senior  years  at 
Vanderbilt,  with  plans  to  continue  their  education  in  graduate  school.  Other 
personnel  affiliated  with  the  project  were  Patricia  C.  Tell inghuisen  as 
Research  Associate  (part  time)  and  myself. 
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Diatomic  rotation -vibration  partition  functions  ate  evaluated  through  a  technique  which  utilizes  semiclassical  phase 
integrals  to  estimate  the  eigenvalues.  For  the  assessment  of  the  contributions  from  the  bound  and  metastable  levels,  the 
method  is  computationally  efficient  and  virtually  exact  for  any  potential  at  any  temperature. 


1.  Introduction  the  spectroscopic  constants.  The  latter  problem  is 

particularly  severe  near  the  dissociation  limit,  where 

In  statistical  mechanics  the  traditional  starting  point  the  centrifugal  distortion  constants  are  important,  and 

for  evaluating  the  internal  partition  function  qur  of  a  where  it  is  also  necessary  to  verify  that  the  calculated 

diatomic  molecule  is  the  combination  of  the  rigid-rotor  levels  actually  exist  in  the  discrete  spectrum.  Moreover, 

model  for  rotation  and  the  harmonic  oscillator  for  vi-  when  discrete  levels  near  dissociation  are  important,  so 

bration  [I],  The  rigid -rotor/harmonic  oscillator  (RR /  are  metastable  levels  (rotational  resonances)  and  even 

HO)  partition  function  may  then  be  corrected  for  ef-  the  vibrational  continuum,  as  Mies  and  Julienne  have 

fects  of  anharmonicity  and  vibration-rotation  interac-  noted  in  an  elegant  treatment  of  the  diatomic  mole- 

tion  through  approximations  [2]  which  are  quite  cuie  in  equilibrium  with  its  parent  atoms  [3], 

reliable  for  temperatures  low  enough  that  only  the  Diatomic  spectroscopic  constants  are  customarily 

first  few  vibrational  levels  contribute  significantly  to  expressed  as  polynomials  in  u  +  j  (4j,  which  is  a  theo- 

qvr.  However,  as  kT  becomes  large  compared  to  the  retically  appropriate  representation  based  on  an  expan- 

vibrational  energy,  hccoe,  it  becomes  necessary  to  in-  sion  of  the  potential  energy  curve  about  its  minimum 
elude  levels  beyond  the  range  of  validity  of  these  ap-  [5].  Near  the  dissociation  limit  these  polynomials  are 

proximations.  In  principle  one  can  always  calculate  not  well  behaved,  exacerbating  the  abov^mentioned 

the  “exact”  qvr  from  an  explicit  sum  of  Boltzmann  problems  of  equating  qvr.  Recent  work  has  shown 

factors  over  all  energy  levels,  using  spectroscopic  con-  that  near-dissociation  expansions  —  which  contain  as 
stants  to  evaluate  the  latter.  As  long  as  only  low-u  argument  uD  —  u,  where  uD  is  the  non-integer  vibration- 

levels  are  important,  this  approach  is  again  satisfactory  al  quantum  number  at  dissociation  -  provide  a  theo- 
and  straightforward  to  apply,  because  the  partition  retically  correct  approach  to  dissociation,  with  effi- 

function  converges  rapidly  to  its  final  value  with  in-  ciency  comparable  to  that  of  the  traditional  polynomials 

elusion  of  at  most  several  hundred  rotational  levels  for  in  v  +  j  [6-9].  Thus  near-dissociation  expansions  offer 

each  of  the  several  low-lying  v  levels.  However,  with  considerable  advantage  in  the  calculation  of  partition 

increasing  T  it  is  necessary  to  include  a  rapidly  in-  functions  when  levels  near  dissociation  are  important, 

creasing  number  of  states  in  the  sum,  and  certain  prob-  However  one  must  still  devise  procedures  for  deter- 

lems  may  arise.  First,  of  course,  there  is  the  obvious  in-  mining  the  range  of  existing  levels  for  various  values  of 

convenience  of  having  to  sum  over  many  more  levels.  the  rotational  quantum  number/,  and  for  checking  the 

But  more  importantly,  the  range  of  included  levels  validity  of  the  expressions  used  to  evaluate  them, 

may  exceed  the  range  of  the  spectroscopic  observations,  In  an  alternative  approach  one  may  start  with  a  po- 
so  that  the  levels  are  no  longer  correctly  represented  by  tential  curve  UU{R)  for  the  molecule  in  question  and 
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numerically  evaluate  the  eigenvalues  of  the  effective 
potential. 

{/yW=W  +  W+l)/^t  (1) 

where  (1  is  a  constant  containing  the  inverse  of  the 
reduced  mass*/  [4  J.  The  rotationless  potential  U0(R) 
can  be  calculated  from  the  spectroscopic  constants 
by  the  RK.R  method,  which  as  generally  applied  is  a 
first-order  semiclassical  method,  but  which  nonetheless 
yields  potentials  which  usually  are  valid  quantum 
mechanically  to  =1  cm-1  [10].  Furthermore  the  po¬ 
tential  Uq(R)  so  determined  can  be  extrapolated  to  the 
dissociation  limit  Oe  as  reliably  and  more  easily  than 
can  the  various  spectroscopic  constants.  And  the  infor¬ 
mation  about  the  existence  of  levels  near  and  above 
De  is  directly  available  from  the  potential:  For  a  given 
/,  energies  above  the  rotational  barrier  (Rh,  £h)  belong 
to  the  continuum,  while  energies  below  the  barrier  are 
discrete  (E  <De)  or  metastable  (E>De). 

Methods  for  calculating  the  eigenvalues  and  wave- 
functions  of  the  effective  potential  Uj(R )  have  been 
available  for  decades  [11—14]  and  can  be  accurate  to 
better  than  I0~3  cm-1 .  However,  for  a  heavy  molecule 
the  numerical  solution  of  the  Schrodinger  equation  for 
all  levels  can  be  quite  time-consuming,  and  also  quite 
tedious  to  apply  near  Z?c,  where  the  potential  flares 
out  gently  to  approach  the  asymptote.  Metastable 
levels  require  special  considerations  [3]  but  can  usually 
be  taken  as  discrete  and  evaluated  by  taking  the  outer 
limit  of  integration  as  Rb,  the  location  of  the  barrier 
[15]. 

In  the  present  paper  I  describe  a  method  of  calcu¬ 
lating  partition  functions  directly  from  potential  curves, 
in  which,  however,  the  eigenvalues  are  obtained  from 
the  semiclassical  phase  integral, 

*2 

+  =  ( 8p)W  /  [E  -  Uj(R)\V'  <1R  .  (2) 

This  equation  yields  directly  the  (non-integral)  vibra¬ 
tional  quantum  number  v  for  any  specified  E  <  ffb, 
through  the  indicated  integral  of  the  radial  momentum 
between  the  classical  turning  points, /?j  and/?2-The 
advantages  of  this  method  are  several:  (I)  It  is  orders 
of  magnitude  faster  than  numerical  solution  of  the 
Schrodinger  equation,  since  the  integrals  in  eq.  (2) 
can  be  evaluated  with  part-per-miUion  accuracy  using 
as  few  as  16  points  [  16] ;  (2)  it  is  straightforward  to 
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apply,  since  the  range  of  integration  is  a  precisely  de¬ 
fined  function  of  the  energy;  (3)  when  used  with  an 
RKR-produced  U$(R),  it  is  “exact”,  since  it  is  simply 
the  inverse  of  the  first-order  semiclassical  RKR  method 
[16].  ,4s  an  illustration  of  these  points  1  note  that  in 
one  of  the  test  cases  described  below,  I  evaluated  more 
than  50000  eigenvalues  and  calculated  partition  func¬ 
tions  at  five  temperatures  -  using  double  precision 
arithmetic  and  less  than  four  minutes  of  CPU  time  on 
a  DEC-1099  computer.  Thus  this  method  is  clearly 
a  practicable  one.  Moreover,  although  the  present  cal¬ 
culations  deal  with  just  the  discrete  and  metastable 
contributions  to  qur,  very  similar  computational  meth¬ 
ods  can  be  used  to  evalute  the  continuum  contributions 
[3]- 

2.  Computational  method 

The  main  trick  which  makes  the  present  semiclas¬ 
sical  approach  feasible  is  the  use  of  the  very  efficient 
Gauss-Mehler  quadrature  to  evaluate  the  phase  inte¬ 
gral  in  eq.  (2)  [16,17]: 

R2  j 

[2/C^2  -  ^l)I  /  [E-Uj(R)\'VdR~  j' F(x  )  dx 

*i  -1 

1  n 

=  /(I  -x-)U2Fw(x)  dx  *  2 HiFJxJ  .  (3) 

-1  1=1 

Note  that  for  a  purely  harmonic  potential,  Fw(x )  is  a 
constant,  and  this  quadrature  becomes  exact  for  any  n. 
The  weights  {//,-}  and  abscissa  {jc;-}  are  simple  tr:gono- 
metric  functions  [18],  which  for  efficiency  are  calculat¬ 
ed  once  at  the  outset  and  stored  for  further  use.  To  il¬ 
lustrate  the  power  of  this  quadrature,  I  show  in  tables 
l  and  2  the  vaiuates  of  t^,  the  vibrational  quantum 
number  at  the  energy  of  the  barrier  fb,  for  selected/ 
levels  of  the  potentials  used  to  approximate  I,  and 
Ar2  in  the  test  calculations  discussed  further  below. 

It  is  clear  that  a  16-point  quadrature  is  more  than 
adequate  in  every  case;  in  fact  most  of  the  time  even 
four  points  are  sufficient!  In  this  regard  it  is  wonh 
noting  that  the  evaluation  of  is  probably  the  most 
severe  test  of  the  quadrature,  since  the  potential  is 
much  farther  from  harmonic  over  the  R  range  needed 
for  this  calculation  than  it  is  for  levels  near  the  ir.ini- 


CHEMICAL  PHYSICS  LETTERS 


5 


]  0 


Volume  102,  number  1 


CHEMICAL  PHYSICS  LETTERS 


1 1  November  1 983 


Table  1 

Numerical  evaluation  of  the  vibrational  quantum  number  at  the  enettty  ot  the  rotational  barrier,  for  the  Morse  potential  approxima¬ 
tive  I; 


J 

vb 

*l 

*b 

Eb 

n  -  4 

n  =  8 

vO 

M 

C 

n  =  32 

n  =  64 

ton 

109.38155 

109.43360 

109.43296 

109.43291 

109.43290 

2.29872 

6.96471 

376.78 

12546.76 

400 

69.50184 

69.49030 

69.48912 

69.48903 

69.48902 

2.37733 

4.89839 

5869.08 

13882.34 

700 

18.39314 

18.38709 

18.38651 

18.38647 

18.38646 

2.64414 

3.85477 

17042.85 

18670.06 

Table  2 

Numerical  evaluation  of  the  vibrational  quantum  number  at  the  energy  of  the  rotational  barrier,  for  the  6-12  potential  approxima' 
ting  Ar  2 


J 

Vb 

ft  -  4 

n  »  8 

n  -  16 

n  =  32 

*1 

Rb  E 

m 

Eb 

1 

7.14407 

9.3S068 

9.48636 

9.48640 

3.34980 

31.66656 

0.12 

100.01 

10 

7.61419 

7.70673 

7.70676 

7.70676 

3.36076 

11.62462 

6.55 

100.47 

30 

4.05543 

4.05558 

4.05551 

4.05551 

3.45029 

6.76853 

54.55 

111.33 

mum.  Also,  near  £  =  £b  a  given  error  in  u  may  repre¬ 
sent  a  relatively  small  error  in  the  estimated  eigenval¬ 
ues,  which  anyway  to  not  need  to  be  evaluated  with 
spectroscopic  accuracy  for  the  calculation  of  qvr.  For 
example,  the  two  mgnest  levels  in  the  effective  poten¬ 
tial  for/  =  1  in  table  2  are  less  than  1  cm-1  apart. 

A  second  procedure  I  have  used  to  save  computa¬ 
tional  time  takes  advantage  of  the  fact  that  eq.  (2) 
gives  v  as  a  continuous  function  of  E.  Thus  one  can 
set  up  a  grid  of  (u,  £ )  values  spanning  the  energy  region 
£m  <  £  <  £b ,  where  £m  is  the  energy  at  the  minimum 
of  Uj(R),  and  then  simply  interpolate  over  this  grid  to 
determine  the  eigenvalues.  In  the  test  calculations  dis¬ 
cussed  below,  as  few  as  seven  points  over  a  grid  linear 
in  £  servedTo  determine  the  energies  of  the  important 
low-u  levels  within  0.01  cm-1 !  (It  is  important  to  in¬ 
clude  the  two  end  points,  (um  =  -0.5, £m)  and  (tib, 
£b)  in  this  grid.)  With  so  few  points  the  levels  in  the 
highest  £  segment  were  occasionally  in  error  by  =10 
cm*1 ;  however  at  low  T  these  levels  are  unimportant, 
and  at  high  T  the  emir  is  insignificant  in  its  effect  on 
qvr.  (At  high  T  it  is  the  number  of  bound  levels  and 
their  degeneracies  that  are  important,  not  their  ener¬ 
gies.) 

With  the  above  points  in  mind,  the  computational 


algorithm  may  be  summarized  as  follows:  (1)  For  a 
given/,  construct  the  effective  potential  Uj(R),  and 
locate  the  minimum  (£m,  £m)  and  barrier  (/?b,  £b); 
(2)  subdivide  the  energy  range  into  L  equal  intervals 
separated  by  A£  =  (£b  -  £m)/£,  and  evaluate  v  at 
the  £  +  1  resulting  values  of  £ ;  (3 )  calculate  the  eigen¬ 
values  £u  for  v  -  0,  1,  2, ...  by  interpolation,  and  eval¬ 
uate  their  contributions,  (2/  +  1)  exp  ( -EJkTi ),  to 
qvr  at  each  of  several  specified  temperatures,  7);  (4) 
repeat  the  procedure  for/  =  0,  1 , 2, ...  until  there  is 
no  bound  level  in  the  effective  potential,  which  occurs 
at  or  slightly  below  the  /  value  wiuch  gives  no  mini¬ 
mum  in  UjiR).  This  procedure  exhausts  the  discrete 
and  metastable  contributions  systematically,  so  that 
only  the  continuum  contributions,  which  must  be  eval¬ 
uated  by  other  means  [3],  remain.  Note  that  the  dis¬ 
crete  and  metastable  contributions  can  be  accumulat¬ 
ed  separately  in  step  (3)  if  that  information  is  of  in¬ 
terest. 


3.  Test  calculations 

I  have  tested  this  method  on  three  potentials  —  a 
Morse  approximation  to  the  ground  state  of  h  (u>e 
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=  ; !  4  cm-1,  £>e  =  1:500cm-1,  Rc  =  2.667  A),  a 
Morse  curve  for  a  typical  “light"  molecule  (q  =  8  umu. 
D.,  =  50000  cm”1.  coe  =  :00Ocm->,jRt,  =  1.2  A),  and 
a  6-1 2  approximation  for  Am  (/a  =  20  amu.  C6 
=  -5.65  X  105  A6  cm”1,  C!2'=  7.98  X  108  A12  cm'1, 
for  which  Re  =  3.76  A  and  De  =  100.01  cm-1 ).  In 
each  case  I  evaluated  the  phase  integral  in  (2)  over 
4,  8,  16, ...  points  until  consecutive  estimates  agreed 
within  An  =  10-4.  For  the  interpolation  grid  L-  6 
proved  adequate;  and  4-point  and  6-point  ^agrangian 
interpolations  gave  identical  values  of  q ur.  As  noted 
earlier,  the  most  computationally  demanding  of  these 
tests  was  the  I->  calculation,  in  which  some  5  X  104 
eigenvalues  were  estimated  over  the  range /  =  0—816. 

In  this  case  I  experimented  with  one  other  nme- 
savmg  measure  -  summing  over  /  =  0,  \f,  ZM,  3 M, 
etc.,  with  multiplication  of  the  Final  sum  by  M.  (The 
motivation  for  this  procedure  is  the  observation  that 
for  homonuclear  molecules  having  zero  nuclear  spin, 
alternate  levels  are  rigorously  absent  from  the  sum.) 

The  results  of  the  calculations  for  the  Morse  poten¬ 
tial  are  compared  with  RR  HO  and  rigid-rotor/Morse 
oscillator  (RR/MO)  values  as  functions  of  Tin  Fig.  1. 
(Note  that  the  latter  is  not  the  same  as  the  full  rota¬ 
tion-vibration  calculation  for  the  Morse  potentials  in 
my  tests.)  As  expected  the  RR/HO  results  are  poorer 
than  the  RR/MO  values,  since  the  latter  correctly 


Eli!.  1 .  Percent  error  in  various  estimates  of  </(ir  for  Ij  fa,  b,  c) 
an  J  Oj  it,  d )  as  functions  of  T.  Curves  (a)  and  (c)  are  for  the 
RR  HO  model,  curves  (b)  and  (d)  for  the  RR/MO  model. 
Curve  (e(  shows  the  percent  error  (x  100)  in  the  "exact”  cal¬ 
culation  fur  M  -  5. 


200  600 
T1K) 


Fig.  2.  Comparison  of  partition  functions  for  Am  calcu'ated 
using  the  RR/HO  model  [curve  (a)],  the  RR/MO  model  lb), 
and  the  “exact”  result  (c).  .Also  shown  is  the  metastarde  con¬ 
tribution  to  the  last  of  these  (d). 

take  into  account  the  vibrational  anharmonicitv  Both 
models  fall  increasingly  below  the  exact  results  as  T 
increases;  this  is  attributable  to  the  omission  of  the  vi¬ 
bration-rotation  interaction,  since  qr  for  the  R  R  mod¬ 
el  is  proportional  to  the  inverse  of  the  rotational  con¬ 
stant  Bv,  and  the  latter  decreases  with  increasing  u  The 
reason  for  the  anomalously  low  values  at  low  T  tot  the 
“02"  tests  are  ( 1 )  the  use  of  Bc  instead  of  j50  in  qr , 
and  (2)  a  5  cm-1  error  in  the  zero-point  energy  ;n  the 
case  of  the  RR/HO  calculation.  With  these  changes  the 
results  for  both  molecules  would  be  much  closer  tc  ex¬ 
act  at  /=  200  K. 

At  temperatures  much  higher  than  those  cove  rut;  in 
Fig.  1,  both  the  RR/HO  and  RR/MO  results  must  even¬ 
tually  exceed  the  exact  qvr.  This  is  evident  in  fig.  2, 
where  the  comparison  is  given  for  the  A m  test  case. 

At  high  T  the  exact  value  approaches  a  constant,  in 
recognition  of  the  existence  of  a  Finite  number  of 
bound  and  metastable  levels.  Both  the  RR/HO  -rid 
RR/MO  values  continue  to  increase  with  T  due  to  the 
direct  proportionality  on  Tin  qr  for  the  RR  modci, 
however,  the  RR/MO  model  correctly  includes  tire 
Finite  vibrational  spectrum,  so  is  somewhat  better  at 
high  T.  Of  course  it  is  precisely  the  region  where  the 
“exact”  qvr  is  leveling  off  that  the  continuum  contri¬ 
butions  become  important  (3],  so  the  tatter  should 
really  be  included  for  a  proper  comparison  (which. 
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however,  is  beyond  the  scope  of  this  work). 

As  shown  by  curve  (e)  in  fig.  1 ,  the  inclusion  of 
only  every  fifth /in  the  rotational  sum  for  I2  leads  to 
an  error  which  is  probably  negligible  everywhere,  ex¬ 
cept  possibly  at  low  T.  Use  of  M  =  2  gave  about  a  fac¬ 
tor  of  four  smaller  error.  The  error,  however,  is  syste¬ 
matic,  and  it  is  clearly  a  function  of  the  particular 
molecule  and  temperature  in  question,  being  smaller 
for  large  q  and  T.  In  any  event  it  appears  that  the 
calculations  are  already  efficient  enough  to  permit 
evaluation  of  the  full  rotational  sum  in  most  cases, 
without  excessive  use  of  CPU  time. 


4.  Conclusion 

The  Morse  oscillator  is  a  workhorse  in  diatomic 
calculations  requiring  a  “realistic"  potential.  Indeed  it 
has  been  used  previously  in  test  calculations  of  parti¬ 
tion  functions  (see,  e.g.  refs.  [19—22]).  However, 
many  such  studies  have  treated  just  the  vibrational 
partition  function  rather  than  the  full  rotation-vibra¬ 
tion  partition  function,  when  in  fact  it  is  the  latter 
that  is  needed  for  thermodynamic  applications.  The 
scheme  I  have  described  is  applicable  to  any  potential 
at  any  temperature.  Of  course  it  does  require  a  poten¬ 
tial  and  will  obviously  be  limited  by  uncertainties  in 
the  potential.  And  it  is  of  limited  use  in  cases  such  as 
that  of  Ari  treated  here,  where  the  continuum  contri¬ 
butions  to  qvr  are  needed  to  account  for  the  total 
“AT2”  (which  may  be  negative!)  in  equilibrium  with 
Ar  atoms.  However,  even  in  cases  such  as  this,  it  can 
be  used  to  assess  total  concentrations  of  bound  and 
metastable  molecules,  or  concentrations  in  selected 
bound  or  metastable  (u,  /)  levels. 

Originally  I  developed  the  present  semiclassical 
method  for  the  purpose  of  testing  partition  functions 
calculated  by  the  classical  method  [3,23],  However, 
in  the  course  of  conducting  the  present  calculations, 
it  has  become  clear  to  me  that  the  method  has  value 
in  its  own  right.  Since  most  of  the  computational  time 
is  invest^wed  in  the  determination  of  the  eigenvalues, 
one  can  evaluate  qvr  at  a  large  number  of  temperatures 
in  a  single  pass,  at  little  additional  expense.  (Sums  in¬ 
cluding  E  and  E *  under  the  summation  sign  can  also 
be  evaluated  directly  rather  than  estimated  from  tem¬ 
perature  derivatives.)  For  the  classical  partition  func¬ 
tion,  on  the  other  hand,  a  completely  new  (albeit 


simple  and  fast)  calculation  is  required  at  each  T. 

To  the  best  of  my  knowledge,  the  simple  method 
I  have  described  here  has  not  previously  been  described 
in  the  literature  and  is  not  known  or  used  by  thermo- 
dynamicists.  However,  since  completing  the  present 
calculations,  I  have  discovered  that  a  very  similar  proce¬ 
dure  was  used  35  years  ago  by  Woolley  et  al.  [24],  to 
estimate  contributions  from  high  (u,/)  levels  to  parti¬ 
tion  functions  of  H2  at  high  temperatures.  These 
authors  even  used  an  RKR-like  procedure  to  determine 
the  width  of  their  rotationless  potential  as  a  function 
of  v.  However,  their  work  was  completed  before  the 
advent  of  modem  computers,  and  it  appears  the  authors 
did  not  recognize  either  the  great  accuracy  of  the 
semiclassical  phase  integral  or  the  efficiency  with  which 
it  can  be  evaluated  using  the  Gauss— Mehler  quadrature. 
Accordingly  their  method  seems  to  have  been  ignored 
in  subsequent  work.  Now  it  is  clear  that  techniques 
based  on  the  semiclassical  phase  integral  can  be  useful 
in  any  calculation  where  levels  near  dissociation  are 
important,  since  they  can  yield  directly  the  number  of 
bound  and  metastable  levels  and  their  energy  range 
(tables  1  and  2),  as  well  as  their  contributions  to  qvr 
(fig-  2). 
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The  B  * X emission  spectrum  is  analyzed  for  the  single  isotopomer  1 ,6Xe  WF.  Thirteen  v'-v"  bands 
spanning  v'  =•  0-4  and  v"  =  0-5  have  been  rotationally  analyzed.  The  analysis  yields  greatly 
improved  spectroscopic  constants,  permitting  precise  identification  of  most  of  the  features  in  the 
XeF  laser  spectrum.  No  evidence  is  found  for  C— *X  emission  in  the  normal  sense.  However,  it  is 
possible  that  several  previously  unassigned  features,  including  one  or  more  reported  laser  lines, 
constitute  C  -X  emission  enabled  through  intensity  borrowing  from  B-X  by  B-C  perturbative 
mining. 

PACS  numbers:  33.20. Kf.  33.70.  -  w.  35.20. Pa.  35.20.Sd 


The  diatomic  rare  gas  halides  burst  into  prominence 
eight  years  ago  as  lasants  in  a  new  class  of  excimer  laser.  At 
that  time  these  molecules  were  virtually  unknown  spectro¬ 
scopically.  Their  essential  properties  were  rapidly  deter¬ 
mined  from  theoretical  and  experimental  studies.'  In  fact, 
most  of  the  lasing  B  -X  transitions  turned  out  to  be  spectro¬ 
scopically  "simple,"  the  only  structure  in  their  bound-free 
emission  spectra  being  Franck-Condon  structure.  Only  in 
the  cases  of  XeF.  XeCI,  and  XeBr  has  discrete  structure  been 
observed.  (In  the  last  case  the  discrete  structure  is  just  barely 
discernible  and  has  not  yet  been  analyzed.2) 

The  most  interesting  RgT  molecule  from  the  spectro¬ 
scopic  standpoint  land  one  of  the  more  interesting  high-pow¬ 
er  laser  candidates)  is  XeF.  The  B-*X  and  D  +X  transitions 
in  XeF  were  analyzed  in  emission '  '  and  absorption6  shortly 
after  the  XeF  laser  was  reported.  Nozzle-beam  studies7  cor¬ 
roborated  the  sketchy  rotational  analysis  of  B-X  in  the  emis¬ 
sion  work  and  gave  additional  information  about  high  v'  and 
low  v"  levels.  Very  recently  the  third  ion-pair  excited  state, 

C 1 3/2 1,  has  been  characterized  through  the  weak  C+—X  tran¬ 
sition/'1  However,  the  spectroscopic  description  of  the  B-X 
transition  remains  marginal  for  the  identification  of  the  las¬ 
ing  transitions,  mainly  because  of  ambiguities  stemming 
from  isotope  effects  in  the  spectrum  of  "natural"  XeF,  which 
contains  predominantly  i:yXe  (26%),  IMXe  (21%),  and 
1  ,;Xe(27%).  To  remedy  this  situation  we  have  been  studying 
the  emission  spectrum  of  the  single  XeF  isotopomer: 

TABLE  I  Spectroscopic  constants  from  band-by-band  rotational  analysis  of  B 
defined  in  Ref  5 


i  tr>xe  i«p  worjc  js  now  sufficiently  completed  to  permit 
reliable  assignment  of  most  of  the  features  in  the  laser  spec¬ 
trum. 

The  emission  spectrum  was  photographed  in  the  region 
3300-3600  A,  using  equipment  and  procedures  similar  to 
those  described  previously. The  Tesla  discharge  sources 
were  charged  with  1-1.5  Torr  SFh,  2-3  Torr  ,  ,6Xe  (Mound 
Laboratory,  95.2%  isotopic  purity),  and  140-170  Torr  Ar. 
(These  pressures  of  SFh  and  Xe,  with  400-600  Torr  Ar,  were 
found  to  give  optimum  intensity;  however,  the  lower  Ar 
pressure  was  employed  to  keep  line  broadening  negligible.) 
The  photographic  plates  had  a  reciprocal  dispersion  of  1.35 
A/mm  and  a  resolving  power  greater  than  200  000.  The 
spectra  were  calibrated  with  Fe  atomic  lines.  We  estimate 
that  sharp,  unblended  lines  are  measured  with  an  absolute 
accuracy  of  +  0.02  cm  which  is  comparable  to  the  stan¬ 
dard  deviations  obtained  in  our  least-squares  fits,  described 
below. 

To  date  we  have  rotationally  analyzed  13  v'-v"  bands 
spanning  v'  =  0-4  and  v"  =  0-5.  In  some  cases  (particularly 
the  strong  bands  from  v'  =  0)  we  are  able  to  extend  the  as¬ 
signments  to  N "  >  70.  However,  when  the  data  are  fitted  to 
the  Hamiltonians  given  for  the  B  and  X  states  in  Ref.  4,  the 
quality  of  the  fit  deteriorates  rapidly  for  N  >  45.  This  prob¬ 
lem  was  noted  earlier  and  is  still  under  investigation.  In  the 
meantime,  to  obtain  optimal  constants  for  the  low-A  lines, 
we  have  fitted  only  the  assignments  for  N  <  41 .  Results  of  the 

-X  spectrum  of 1  “’Xe  'T.  All  quantities  are  given  in  cm  Constants  arc  as 
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'Analyzed  in  Ref  5  Present  results  indicate  a  probable  calibration  error  of  081cm  '  in  the  data  in  Table  I  of  Ref  5 
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hand-by-band  tits  arc  summarized  in  Table  I.  in  which  the 
constants  arc  as  defined  in  Refs.  4  and  5.  The  constants  for 
each  band  represent  typically  20-25  assigned  lines  in  each  of 
the  four  branches  (P,  ,  Pr  R  ,  Rt ),  although  the  assignments 
are  somewhat  sparser  for  those  bands  involving  r"  =  4  and 
5,  in  which  the  rotational  structure  is  very  congested. 

To  obtain  the  usual  spectroscopic  expansion  coeffi¬ 
cients  for  the  fi-.Vband  system,  we  proceeded  as  follows.  The 
estimated  band  origins  from  Ref.  5  were  augmented  by  the 
band  origins  from  Table  I  (heavily  weighted  in  accord  with 
their  much  greater  precision)  and  were  fitted  directly  to  the 
usual  double  polynomials  in  [p\v’  +  1/2|]  and  [p\ v"  +  1/ 
2)],  where  p  is  the  isotopic  factor,  related  to  the  reduced 
masses  by  p  =  (//,/«)'  We  take  1  w‘  XeF  to  be  the  reference 
molecule  (fir  I.  The  data  from  Ref.  5  were  taken  to  represent 
the  hypothetical  “average"  molecule,"1  'XeF,  for  which 
p  =  1.00215.  | Attempts  to  give  the  various  features  more 
specific  isotopic  assignments  actually  gave  poorer  variances, 
hence  were  abandoned.)  To  fit  these  data  adequately,  it  was 
necessary  to  use  seven  vibrational  parameters  for  the  X  state 
|cf.  five  in  Ref.  5).  This  is  indicative  of  the  anomalous  nature 
of  the  XeF  ground  state,  since  the  number  of  vibrational 
parameters  is  only  three  less  than  the  number  of  levels  repre¬ 
sented  (0-9).  For  the  B  state  three  parameters  gave  a  slightly 
better  fit  than  two,  with  the  data  spanning  v'  =  0-15. 

To  obtain  the  rotational  and  spin-splitting  constants, 
w  e  employed  a  correlated  fit  of  the  relevant  results  in  Table  I 
and  their  associated  variance-covariance  matrices,  using 
methods  like  those  outlined  by  Albritton  et  al."  and  used 
previously  by  us."  Flere,  again  we  found  that  the  A-state 
rotational  constants  were  not  well  overdetermined,  as  five 
parameters  were  needed  to  represent  the  six  observed  levels 
0-5.  iThe  five  expansion  parameters,  however,  were  better 
than  a  direct  fit  to  the  six  individual  B,''  values,  and  they 
gave  virtually  identical  results.)  Two  parameters  were  ade¬ 
quate  for  B,  and  three  for  D,  ".while  three  were  required  for 
both  SB, '  and  aB,  "  [also  expressed  as  polynomials  in 
It’’  +  1/2)  and  It'"  4-  1/2)].  The  “extra”  centrifugal  distor¬ 
tion  parameter  b,  which  was  introduced  in  Ref.  5,  showed 
stronger  correlation  with  v'  than  with  e",  so  was  expressed  as 
a  linear  function  of  In'  +  1/2).  These  results  are  summarized 
together  with  the  vibrational  parameters  in  Table  II.  It  is 
worth  noting  that  our  results  are  similar  to,  but  not  com¬ 
pletely  consistent  with,  the  only  other  rotational  analysis  of 
this  system,  that  of  Ref.  7.  For  example,  the  estimate  of  R,  " 
in  the  latter  work  was  2.321(2)  A. 

Although  the  constants  in  Table  II  are  strictly  valid  for 
only  the  '  ‘'’XeF  isotopomer,  we  expect  them  to  be  applicable 
to  the  other  isotopomers  as  well,  at  least  to  a  good  approxi¬ 
mation.  To  use  them  in  this  way,  one  replaces  (c  +  1/2)  by 
plr  f  1/2)  everywhere,  and  k  [^J{J  +  1)  or  A  (A  +  1),  see 
Refs.  4  and  5]  by /rtf,  where/?  =  1.00332  for  ':<JXeF,  1.00234 
for  "'XeF,  1.00186  for  "’XeF,  and  1.00092  for  '"XeF.  We 
assume  that  the  spin-splitting  parameters  8  and  a  show  no 
specific  isotope  dependence.  Since  these  parameters  are  giv¬ 
en  in  Table  II  in  the  coupled  forms  SB, '  and  aB,  " ,  in  which 
they  were  determined  in  the  fits,  we  have  extracted  and  pre¬ 
sented  in  Table  III  theexplicit  values.  It  should  be  noted  that 
these  constants  are  much  less  precisely  determined  than  are 
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'The  given  quantities  are  coefficients  of  polynomials  in  ii  *  l/2t.  with  nu¬ 
merical  subscripts  indicating  the  power  of  the  argument  Where  given,  tilt- 
figures  in  parentheses  represent  standard  errors,  in  terms  of  final  digits 
b  Vibrational  parameters  valid  for  r  =  0-15,  r"  0-4;  rotational  param¬ 
eters  valid  for  r'  -  0-4,  r"  =  0-5  The  parameters  are  given  to  sufficient 
precision  to  permit  recalculation  of  original  results  within  experimental  er¬ 
ror 

•  Quantities  not  given  le  g..  D. '  and  //,  "|  remain  as  estimated  in  Table  V  of 
Ref.  5 

'Coefficients  of  polynomials  for  aB,  "  and  i SB, ' .  respectively 

B,  ’  and  B,  ",  In  fact,  8  and  a  show  very  high  negative  correla¬ 
tion,  so  that  the  quantity  which  is  well  determined  is  the  sum 

SB,  '  +  aB,  " . 

The  strongest  features  in  the  XeF  laser  spectrum  are 
groups  of  lines  near  3511  and  3532  A,4"  with  the  former 
predominating  under  some  circumstances  and  the  latter  un¬ 
der  other.1'’  IK  The  present  results  solidify  the  earlier,  more 
tentative  assignment  of  these  features.4  For  example,  we 
show  in  Fig.  1  the  high-resolution  spectrum  near  3511  A. 
with  calculated  rotational  assignments  for  the  1-4  and  0-2 
bands.  Even  under  the  present  very  high  resolution,  the  rota¬ 
tional  structure  in  this  region  is  only  marginally  resolved.  In 
a  high-pressure  laser  operating  on  “natural"  XeF,  isotope 
shifts  and  pressure  broadening  produce  further  blurring  of 
this  structure.4 


TABLE  III.  Spin-splitting  parameters  for  .V  and  B  states  of  XeE 
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FIG  I  Densitometer  tracing  of  the  emission  spectrum  of 1  "XeF  near  3511 
A,  showing  rotational  branch  structure  in  the  1-4  band  Also  shown  Ito 
righti  is  the  R.  head  in  the  0-2  band 


In  their  recent  analysis  of  the  CV  X  spectrum.  Helm  el 
a!."'  have  located  the  C  state  near  —  2.45  A  and  28023 
cm  or  about  0.2  A  to  smaller  R  and  789  cm  'below  the  B 
state.  They  have  also  suggested  that  for  thermally  equilibrat¬ 
ed  B  and  C  states,  the  C  * X  emission  could  account  for  as 
much  as  35*T  of  the  total  UV  emission.  We  have  looked  for 
C  -X  bands  in  our  spectra  and  have  found  no  evidence  of 
them.  In  this  regard,  we  note  that  the  strongest  C  - X  bands 
should  be  0-1  and  0-2.  which  art  calculated  to  lie  near 
27  883  and  27  7(X)  cm  '.  respectively.  This  is  in  the  red 
wings  of  the  B  -X  system,  w  here  only  very  weak  emission 
occurs  from  r  0  to  c"  -*  7,  so  C  • X  emission  should  be 
easily  discerned.  Since  our  single  isotope  sources  were  oper¬ 
ated  at  low  pressures,  we  also  reexamined  the  Sandia  spec¬ 
tra.4  w  hich  were  recorded  at  P-  2  atm.  Again  we  found  no 
evidence  of  the  predicted  C  -X  bands.  From  this  result  we 
conclude  that  the  C  »,V  system  must  be  at  least  a  factor  of 
-  30  weaker  than  suggested  in  Refs.  8  and  9. 

The  analysis  of  the  C  state  also  permits  us  to  predict 
perturbations  in  the  B  *.V  spectrum.  For  low  r„,  levels 
V,  r„  *  2  lie  somewhat  below  eH,  catching  up  near 
i  „  3  Also,  because  the  rotational  constants  are  larger  for 

C  than  for  B.  the  rotational  levels  of  C  overhaul  those  of  B 
from  below,  so  that  perturbations  should  occur  for  rota¬ 


tional  quantum  numbers  as  low  as  —  20  in  uu  =  3,  or  as  high 
as  —  70  in  vB  =  0.  So  far  we  have  not  been  able  to  identify 
such  perturbations,  although  the  search  has  not  been  ex¬ 
haustive.  However,  w edo  think  it  likely  that  the  previously 
mentioned  problems  with  the  fit  Hamiltonian  are  due  to  an 
anomalous  centrifugal  distortion  of  the  B  state  due  to  inter¬ 
action  with  the  C  state.  It  is  also  possible  that  some  of  the 
previously  unassigned  features4  may  in  fact  be  C— <-X  emis¬ 
sion  from  C-state  levels  that  borrow  intensity  from  B—*X 
through  B-C  perturbative  mixing.  In  this  regard  it  is  interest¬ 
ing  to  note  that  one  of  the  weaker  laser  lines  reported  by 
Shimauchi  et  a/.15  (at  3502.475  A)  coincides  with  one  of  the 
strong,  unassigned  features  in  Ref.  4.  Thus  it  is  conceivable 
that  the  C  -*X  transition  has  even  been  made  to  lase,  though 
in  a  rather  special  way!  We  are  continuing  work  on  these 
problems. 

This  work  was  supported  by  DARPA  and  ONR. 
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